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ABSTRACT 

We obtained the number counts and the rest-frame _B-band luminosity function of the color- 
selected old passively-evolving galaxies (OPEGs) at z = 1 with very high statistical accuracy 
using a large and homogeneous sample of about 4000 such objects with z' < 25 detected in the 
area of 1.03 deg^ in the Subaru/XMM-Newton Deep Survey (SXDS) field. Our selection criteria 
are defined on the i' — z' and R ~ z' color- magnitude plane so that OPEGs at z = 0.9 — 1.1 
with formation redshift z^=2-10 are properly sampled with minimum contamination by other 
population. The limiting magnitude corresponds to the luminosity of galaxies with Af, -I- 3 at 
z — 0. We made a pilot redshift observations for 99 OPEG candidates with 19 < z' < 22 and 
found that at least 78% (73/93) of the entire sample, or 95% (73/77) of these whose redshifts 
were obtained are indeed lie between z = 0.87 and 1.12 and the most of their spectra show 
the continuum break and strong Ca H and K lines, indicating that these objects are indeed 
dominated by the old stellar populations. The relationship between the observed redshift and 
the color well follows the models used in defining the selection criteria in the consistent manner. 
We found that the surface number density of OPEGs varies 10-30% of the whole average from 
field to field even at a 30-arcmin scale and the poorest field, SXDS-South, has only 65% of the 
richest one, SXDS-East. We then compare our results with the luminosity functions of the color- 
or the morphologically-selected early type galaxies at z = taking the evolutionary factor into 
account and found that the number density of old passive galaxies with magnitude at 

z ^ 1 averaged over the SXDS area is 40-60% of the equivalently red galaxies and 60-85% of the 
morphologically-selected E/SO galaxies at z = depending on their luminosity evolution. It is 
revealed that more than half, but not all, of the present-day early-type galaxies had already been 
formed into quiescent passive galaxies at z = 1. 

Subject headings: cosmology: observations — cosmology; large-scale structure of the universe — galaxies: 
high-redshift — galaxies: photometry 

1. INTRODUCTION 



^ Based on data collected at Subaru Telescope, which 
is operated by the National Astronomical Observatory of 
Japan. 

^Optical and Infrared Astronomy Division, National As- 
tronomical Observatory, Mitaka, Tokyo 181-8588, Japan. 

^European Southern Observatory, Karl-Scwarzschild- 
Str.2, D-85748, Garching, Germany. 

^Subaru Telescope, National Astronomical Observatory 
of Japan, 650 N. A'ohoku Place, Hilo, HI 96720, USA. 

^Okayama Astrophysical Observatory, National Astro- 
nomical Observatory of Japan, Kamogata, Okayama, 719- 
0232, Japan. 

*^ Space Telescope Science Institute, 3700 San Martin 



Galaxies are expected to evolve passively af- 
ter the major star formation ceased by consuming 
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and/or losing a large fraction of their gas. Their 
photometric properties then change only due to 
the aging of the stars inside, which can be traced 
very well by simple stellar evolutionary synthe- 
sis models without suffering large effects by dust 
extinction or any consequence of on-going star- 
formation activity. Only merging or interaction 
with other galaxies may change their properties 
significantly. 

Aragon-Salamanca et al. (1993) first reported 
detection of systematic change of the colors in 
cluster galaxies at 2;=0.5-0.8 which are consistent 
with passive evolution models with formation red- 
shift, Zf, larger than 2. This was followed by a 
large number of works studying passive-evolution 
properties in cluster and field galaxies at interme- 
diate redshift (e.g., Yamada and Arimoto 1995; 
Stanford et al. 1997; Kodama et al. 1998; Zepf 
1997). The expected characteristic colors of pas- 
sive galaxies are also used to search for clusters 
of galaxies at high redshift (e.g., Dickinson et al. 
1995; Yamada et al. 1997; Tanaka et al. 2000; 
Gladders et al. 2003). The survey for extremely- 
red objects (ERO) also sampled a number of old 
passively evolving galaxies (hereafter OPEGs) at 
high redshift (e.g., Cimatti et al. 2000; Daddi et 
al. 2002; Miyazaki et al. 2003; lye et al. 2003). 

Since galaxy colors are easier to be measured 
and modeled than other properties such as size, 
velocity dispersion, spectral line index, or mor- 
phology, it is useful to construct a sample of well- 
defined color-selected galaxies at each epoch, espe- 
cially at high redshift, in order to firmly constrain 
galaxy formation and evolution scenarios. OPEGs 
are indeed suitable population for such a purpose 
since their colors are even more characteristic at 
high redshift and their photometric evolution is 
the simplest. 

In this paper, we focus on the number density of 
OPEGs at high redshifts, since it provides a crit- 
ical test for the build-up of such massive galax- 
ies as predicted by the hierarchical formation of 
galaxies in the CDM Universe (e.g., Kauffmann 
& Chariot 1999). There have been a number of 
attempts in this line measuring the number den- 
sity of high-^ OPEGs and comparing it with lo- 
cal counterparts, such as COMBO-17 (Bell et al. 
2004), K20 (Pozzetti et al. 2003), GDDS (Glaze- 
brook et al. 2004), MUNICS (Saracco et al., 2005), 
and GOODS (Treu et al. 2005). However, the re- 



sults are often inconsistent, and nothing is conclu- 
sive yet, largely due to the limited spatial coverage 
probably. 

To remedy this situation, we here present a re- 
sult from the Subaru/XMM-Newton Deep Survey 
(SXDS, Sekiguchi et al. 2005, in preparation) on 
the number density of the well defined sample of 
OPEGs at 2; = 1. SXDS is designed to have unique 
combination of large area (1.2 deg^) and depth 
[B < 27.5, R < 27, i' < 27, and z' < 26, all 
in AB system), which enables us to construct a 
sample of OPEGs down to a few magnitude be- 
low the characteristic luminosity, L^, even at such 
high redshift. We analyze the distribution of the 
4118 2: ~ 1 OPEG candidates down to z' = 25 se- 
lected from the total of 330,000 SXDS ^'-selected 
galaxies in the SXDS Verl.O catalog (Furusawa et 
al. 2005, in preparation). The large volume cover- 
age of SXDS is crucially important to suppress the 
uncertainties from field-to-field variation. Since 
strong clustering of OPEGs at high redshift is nat- 
urally expected in biased galaxy formation scenar- 
ios, large volume coverage is essential to sample a 
fair fraction of the universe, only with which we 
can discuss their averaged properties of galaxies. 

With the current optical data-set, the highest 
redshift at which we can fully trace the OPEGs 
is 2: ~ 1, since our longest wavelength filter z'- 
band can still sample the light from the wavelength 
longward of the rest-frame 4000A break. Also, at 
lower redshifts, OPEGs tend to be severely con- 
taminated by blue foreground/background popu- 
lations. Therefore, z = 1 is the optimal choice of 
redshift to explore in the current study. We note, 
however, that in the very near future, the SXDS 
field is a subject of deep near-infrared (NIR) imag- 
ing in JHK bands through the UKIDSS consor- 
tium (Warren 2003), which will then enable us to 
trace OPEGs toward much higher redshift up to 
2; ~ 3 with a thorough photometric-redshift anal- 
ysis. 

The model-dependent color selection must be 
calibrated by the detailed spectroscopic observa- 
tions. For the purpose, we also conducted spec- 
troscopic follow-on observations using the Subaru 
telescope. Although it is time consuming and our 
sample is still small, the results show the color se- 
lection criteria works well in the very consistent 
manner. 

The outline of the paper is as follows. In Sec. 2, 
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wc describe our eolor and magnitude criteria to se- 
lect OPEGs and the effects of on-going minor star- 
formation activity, as well as the brief introduction 
to the SXDS data set. In Sec. 3 we describe our 
spectroscopic observations for the selected OPEG 
candidates. We then present the number counts of 
OPEGs and their variation across the SXDS fields 
in Sec. 4. The rest-frame i?-band luminosity func- 
tion is presented and discussed in Sec. 5. For the 
related topics, the color properties of the general 
2 ~ 1 galaxies (i.e., not only OPEGs but also bluer 
galaxies) in SXDS fields are discussed in our previ- 
ous paper, Kodama et al. (2004). The clustering 
properties of the OPEGs will be also discussed in 
the future publication. Throughout the paper, we 
adopt a set of cosmological parameters as follows: 
the matter density 0^=0.3, the cosmological con- 
stant fl\=0.7, and the Hubble constant Hq=7Q 
hro km Mpc~^. 

2. SAMPLE SELECTION 

2.1. Subaru/XMM-Newton Deep Survey 

Subaru/XMM-Newton Deep Survey (SXDS, 
Sekiguchi et al. 2005, in preparation) is a deep 
and wide-field optical/x-ray survey to observe a 
significant volume of the universe from interme- 
diate to high redshift in order to study structure 
formation history of the universe with high statis- 
tical accuracy. The significant part of the area will 
be also observed in various wavelengths from ultra- 
violet to radio. At optical wavelength, the field 
with 1.2 deg2 area centered at (a, 6)={02^18"', 
-05°00") (J2000.0) was surveyed with Subaru 
Telescope (lye et a. 2004) equipped with Suprime 
Cam (Miyazaki et al. 2002) in B, R, i', and z' 
band by late 2003. B- and i?-band magnitudes 
are calibrated to those in Johnson-Cousins sys- 
tem and i'- and 2;'-band filters are to the Sloan 
Digital Sky Survey. The magnitude values are all 
measured in AB system. The survey field consists 
of the five Suprime-Cam field of view, namely, 
SXDS-Center, -North, -South, -East and -West 
(hereafter SXDS-C, -N, -S, -E, and -W fields). 
The SXDS optical source catalog Verl.O has been 
constructed by October, 2004. The detection limit 
of the shallowest part of the SXDS Verl.O data 
set isB = 27.5, R = 27.0, i' = 26.9, z' = 25.8 (5cr 
sky fluctuation with 2" -diameter apertures). The 
seeing sizes of the composite images are 0. 80- 



0."84. The properties of SXDS-^erl.O discussed 
in this paper is summarized in Table 1. The de- 
tail of the data reduction and the photometry will 
be presented in Furusawa et al. (2005, in prepa- 
ration). The colors used in this paper are those 
obtained with 2" -diameter aperture after correct- 
ing the small 0.05 mag) aperture magnitude 
differences between each fields, which is due to the 
slightly different PSF shapes. We do not use the 
data within 300 pix of the edge to avoid possible 
artificial effects. We also avoid the area that is 
affected by the strong CCD blooming effect by 
the bright sources. The final effective area used in 
this paper is 1.025 deg^. 

2.2. Photometric Sample of Old Passively- 
Evolving Galaxies 

We selected a sample of old passively-evolving 
galaxies (OPEGs) in the SXDS field as follows. 
Fig.l and 2 shows that expected colors and mag- 
nitude evolution of OPEGs with 5L* and O.IL* 
(at z = 0) formed at Zf = 10, 5, 3, 2, and 1.5 along 
the observed redshift, which are calculated using 
the Kodama and Arimoto (1997) evolutionary syn- 
thesis models. We then defined the criteria for 
z = 1 OPEGs on the color-magnitude diagram; 
(i) 0.8 <i' - z' < 1.2, and (ii) -0.05z' + 3.01 < 
R- z' < -O.OSz' + 2.49 to pick up OPEGs at 
z = 0.9 — 1.1. The first criteria is to chose red 
galaxies with strong continuum break at around 
8000 A. Note, however, this i' — z' criteria alone is 
not sufficient to isolate the pure passively-evolving 
galaxies. We therefore adopted the second crite- 
rion which brackets the expected color-magnitude 
relation for OPEGs formed at z/ = 2 — 10 and 
observed &t z = 1. Considering the photometric 
errors, we have further 0.05 mag margin at both 
red and blue side of the model loci. 

Fig. 3 shows the z' versus i' — z' color-magnitude 
diagram of the ^'-selected galaxies in SXDS field. 
To avoid crowdedness in the figure, one-third of 
the whole objects, which are randomly selected by 
their positions are plotted. While we do not use 
any morphological criteria to select OPEG candi- 
dates above, the sources with the SExtractor stel- 
larity index smaller than 0.9 are plotted by the 
smaller red dots in Fig. 3 to demonstrate the galaxy 
color distribution. There is a red envelope in i' — z' 
color as expected from spectral evolution of nor- 
mal galaxies. 
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Fig. 4 shows the z' versus R — z' color- 
magnitude diagram for the objects which satisfy 
our criterion (i). The criterion (ii) is represented 
in the figure by the two tilted solid lines. While 
many of those i' — z' color-selected red galaxies are 
distributed in the range defined by the criterion 
(ii), there is notable scatter in their R — z' color. 
Among the 10655 objects with 0.8 < i' - z' < 1.2 
and z' < 25, » 39% (4118) satisfy the criterion 
(ii) while w 25% (2619) have the redder R - z' 
color and « 37% (3918) have the bluer color. We 
selected the objects which satisfy criterion (i) and 
(ii) as the candidates of OPEGs at ^ ~ 1. Their 
sky distribution is shown in Fig. 5. 

The numbers of the detected OPEGs candi- 
dates are summarized in Tabic 2. Complete- 
ness and contamination for these color selection 
of OPEGs are discussed below and the calibration 
by the spectroscopic follow-up observation is pre- 
sented in Section 3. 

2.2.1. Completeness in color selection 

While our criteria arc set so that we can pick 
up OPEGs at ^ = 0.9 - 1.1 with z/ = 2 - 10 com- 
pletely, OPEGs with younger age or galaxies with 
some star-formation activity may run out from the 
color range. We already reported in our previ- 
ous paper, Kodama et al. (2004), how the i' — z' 
and R — z' colors behave along the rcdshift for the 
models with different ratio of old and young (or 
'bulge-to-disk') populations. Prom Fig.l in this 
paper and Fig.l in Kodama ct al. (2004), one 
may see how our criteria work to reject galaxies at 
lower redshift or galaxies with young population. 

Here we consider more subtle cases, namely, the 
effect of a small fraction of young stellar compo- 
nents added to the colors of old passive galaxies. 
For example, if we add the spectrum with on-going 
star formation component to OPEG spectra with 
Zf > 2, R — z' color becomes bluer rapidly. We 
check how much star-formation brings them out 
from our color range for OPEGs. We consider the 
two extreme cases, namely, constant continuous 
star formation and episodic starburst on a pas- 
sively evolving model galaxy with Zf = 5 with 
luminosity at 2; = 0. First, we add the young stel- 
lar population with constant star-formation rate 
started at 2 = 5 to this old population which has 
i' - z' = 1.05 and R - z' = 1.98 at z = 1. We 
found that their R — z' color becomes bluer than 



the limit [R — z' ^ 1.7) when wc add the young 
population with more than 5% of the stellar mass 
while i' — z' color becomes bluer than 0.8 only if we 
add more than 10%. Second, wc also consider the 
cases that we add young 'burst' (namely, episodic) 
stellar population. If we add the burst component 
with age of 0.2 Gyr with mass fraction of 2-3%, 
R — z' becomes bluer than the limit. If the age 
of the burst is as old as 0.5 Gyr, however, even 
adding the 10% mass of the young component docs 
not remove the model out from our criteria. Only 
very young population seen in relatively short time 
scale brings them out from our criteria. 

On the other hand. Fig. 2 shows that old and 
luminous galaxies at only slightly high redshift be- 
comes redder than the limiting R — z color while 
their i' — z' color are still bluer than 1.2. So those 
galaxies with redder R — z' color maybe galax- 
ies with slightly higher than z = 1, which implies 
sharp cut off in redshift to our sample. 

2.2.2. Contamination by dusty red objects 

While previous search for extremely red objects 
(EROs) with R-K >3 (AB) down to if = 22-24 
contains similar population of old passive galaxies 
at around z = 1, it is known that a simple color 
cut also samples dusty red galaxies at intermedi- 
ate or high redshift. For example, Miyazaki et al. 
(2003) reported that in their sample of EROs se- 
lected by the criteria R — Ks > 3.35 (AB) and 
Ks < 22.1 (AB), only 58% are the old galaxies at 
z=l-1.5. With the Calzetti's extinction low, the 
blue objects at 2; ~ 1 — 2 with intrinsic spectra 
/y=constant have colors 0.8 < i' — z' < 1.2 when 
they suffer reddening with 1.0 < E{B — V) < 1.4 
(at z = 1-2). Since objects with E{B -V)> 1.25 
have R—z' > 2.0, some fraction of dusty red galax- 
ies arc thus rejected by our second criteria. Spec- 
troscopic observations may tell more quantitative 
number for the contamination (see Sec. 3) 

We would like to note that we do not use the 
term 'extremely red galaxies' for our sample since 
it is well in the range of the expected galaxy col- 
ors and not 'extreme' in any sense. We did not 
make any star /galaxy separation for the primary 
sample of the OPEG candidates analyzed below, 
either, although we made various tests includ- 
ing/excluding stellar sources selected by various 
threshold of stellarity index in the SXDS catalog. 
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2.2.3. Com,pleteness in observation near the lim- 
iting magnitude 

We evaluated the incompleteness of our galaxy 
sample due to the photometric error along their 
apparent magnitude by using the model galaxies 
as follows. We first created the blank 'sky' im- 
ages after removing the objects detected in the 
SXDS Verl.O catalog. The areas of the objects arc 
substituted by the 'median' sky image constructed 
from three images at the different part of the same 
Suprimc-Cam field. Then wc rondomly put the 
model galaxies with the similar number density 
(of not only OPEGs but the entire 2;'-band se- 
lected objects) as observed on to these 'sky' images 
and counted the number of the objects recovered 
with the magnitude difference less than 0.25 mag 
and their fraction to the input ones. The model 
galaxies have de Vaucouleurs profile, assuming the 
most of the OPEGs are early-type galaxies, and 
their size arc set to roughly follow the Kormendy 
relation. Fig. 6 shows the completeness of such z'- 
selected galaxies in SXDS-C, -N, -S, -E, and -W 
fields and their average. One notes that the eval- 
uated completeness are very similar for the five 
fields; at z' = 25, the completeness of the SXDS-S 
field with the shortest exposure time (11040 sec) is 
67%, that is not far from other fields with 70-73%. 
The detection completeness is higher than 90% at 
z' = 23, which corresponds to the OPEGs with 
« M« -|- 1. We note that the procedure adopted 
here may slightly (a few %) overestimate the de- 
tection completeness since wc did not fully took 
the effects of the source confusion and clustering 
into account. More comprehensive study of the 
detection completeness of the SXDS catalog will 
be discussed in Furusawa et al. (2005, in prepara- 
tion). 

2.2.4- Photometric errors in color 

It is possible that the sample is contaminated 
by the bluer or redder galaxies or becomes incom- 
plete due to the systematic offset of the magni- 
tude zero point as well as the random photomet- 
ric errors at the faintest magnitude range since the 
number of the objects increases rapidly toward the 
fainter magnitude, especially for the blue ones. 

We first evaluate the possible effects of magni- 
tude zero point errors as follows. We shifted the 
color criteria slightly blue or red, and examined 



how many the numbers of the selected galaxies 
change. If we shift the R — z' colors alone 0.05 
mag, either to the red or the blue side, the num- 
ber of OPEG candidates decreases net 4% (14% 
and 15% drop off and 12% and 11% comes in to 
the sample when shifted to the red and the blue 
side, respectively) for the z' < 25 samples, which 
is only twice of the Poisson error of the sample. 
The sample seems stable for the effect. If we shift 
both of the colors i' — z' and R — z' to the same 
direction, however, the number difference becomes 
as large as 14%; if we shift the criteria 0.05 mag 
bluer, the sample increase 10%, and if 0.05 mag 
redder, it decrease 14%. After all, the uncertainty 
of the number densities of the OPEGs evaluated 
by the color selection alone is likely to be 5-15% 
cannot be rejected. This fact reminds us of the 
importance of the spectroscopic calibration for the 
model-dependent color criteria. Our rather pilot 
study, presented in the paper, indicates that the 
relationship between the redshift and i' — z' color 
follows the adopted model tracks very well (see 
Sec. 3 and Fig. 11 below). Larger sample of spec- 
troscopic redshift or the robust photometric red- 
shift for the entire sample after UKIDSS will be 
much helpful in isolating the red-sequence galaxies 
with rather model-independent empirical manner. 

Next, we consider the effects of random photo- 
metric error of the bluer galaxies at the faint end. 
Typical errors in R — z' and i' — z' colors is 0.08 
and 0.2 mag at z = 24.5 and = 25.5, respectively. 
Assuming the Gaussian distribution of the photo- 
metric error with the dispersion of these values, 
we evaluated the possible contamination by cal- 
culating the objects coming in and dropped from 
the criteria. We ignore the effects due to the first 
stage, namely the i' — z' color selection, and con- 
sider the effect on the z' vs. R—z' color-magnitude 
diagram. While the number of the OPEG candi- 
dates at z' = 24-25 and 25-26 are 803 and 743, 
respectively, the net numbers of the contaminants 
are expected to be « and 195, which yields the 
nominal contamination fraction 0% and 26%, re- 
spectively. We here ignore the bias that the prob- 
ability that the red galaxies are dropped out from 
the criteria is slightly larger than that blue galax- 
ies come in due to the larger error in their color 
for the red galaxies. It is seen in Fig. 4 that 
the most of the contaminants are from the objects 
with z' > 25.5. Thus the OPEG sample seems not 
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affected by the effect above 2' = 25, but it can be 
significantly affected by the contamination of the 
bluer objects at the fainter magnitude. 

3. SPECTROSCOPIC SAMPLE 

While we discussed on the possible contamina- 
tion and incompleteness of our OPEG selection 
criteria in the last section, the true selection func- 
tion can be only determined by spectroscopic ob- 
servations. For the purpose, we conducted optical 
spectroscopy at a fraction of the SXDS field. 

3.1. Observations 

In 2003 October and November, SXDS team 
conducted the follow-up spectroscopic observa- 
tions of the optical/x-ray imaging surveys using 
Faint Object Camera And Spectrograph (FOCAS, 
Kashikawa et al. 2002) equipped with Subaru tele- 
scope. We observed more than 20 slit-aperture 
masks during the run. Many mask fields are cho- 
sen so that the interesting x-ray selected sources 
are to be observed but slitlets are also put on the 
OPEG candidates with z' < 23 in such fields. 
If there is enough room, we also put slitlets at 
OPEG candidates with 23 < z' < 24 as well as 
the non-OPEGs with relatively small color differ- 
ence. These target fields are therefore not ho- 
mogeneously selected from the whole SXDS area 
but it is still unbiased for the red galaxies un- 
less there is strong correlation between the x-ray 
sources and OPEGs at z ^ 1. While there are 
in fact some concentration of the x-ray sources at 
z = 1.1 (Akiyama et al. 2005, in preparation), 
most of the other x-ray sources are at different 
redshift. We refer these masks as 'general field' 
masks for the OPEGs. Besides these masks, we 
also observed two masks in order to investigate 
the redshift distribution of the two 'high density' 
regions of red galaxies, CI and C2, discussed in 
Kodama et al. (2004). In these two masks not 
only the red galaxies but also bluer galaxies are 
observed simultaneously. 

The sky distribution of the OPEG candidates, 
for which we observed their spectra, is shown in 
Fig. 7 by the open and filled circles. In total, we 
observed 175 OPEG red galaxies with 19 < z' < 
24, of which 157 are in the 'general field' masks and 
other 18 are in the 'high-density region' masks. 
The 18 sources are located within 3' radius from 
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Fig. 1. — i' — z' and R — z' colors of the model 
galaxies. The passive evolution cases with the for- 
mation epoch Zi?=10, 5, 3, 2, and 1.5 (from red to 
blue) are shown for the galaxies with 5i» (thick 
lines) and O.IL* luminosity at present. 
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Fig. 2. — Same as Fig.l but the colors are plot- 
ted along the apparent magnitude. The filled and 
open squares show the redshift z = 1 and 2. The 
dashed lines are the color criteria for old passively- 
evolving galaxies adopted in this paper. 
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Fig. 3. — z' versus i' — z' color-magnitude diagram 
of the z'-selected galaxies in SXDS field. To avoid 
crowdedness in the figure, one-third of the whole 
objects which are randomly selected by their po- 
sitions are plotted. The red dots are those with 
stellarity index larger than 0.9. Dashed lines indi- 
cate the detection limits in z' and i' bands for the 
shallowest field. The solid lines show the adopted 
color criteria. 



see Fig4.jpg 
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Fig. 4. — z' versus R—z' color-magnitude diagram 
for the objects with 0.8 < i' — z' < 1.2. Dashed 
lines indicate the nominal detection limits in z' 
and i' bands. Solid lines indicate the adopted color 
criteria for the color. 
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Fig. 5. — The sky distribution of the selected 
OPEG candidates with z' < 25 is shown. 
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the center of the regions CI and C2 in Kodama at 
al (2004). In Fig. 7 we also show the locations of 
the six high-density regions in Kodama et al. by 
the circle with 3' radius. 

We used the '150' grism except for the two 
masks we observed with the '300B'grism. Slit 
width is 0.8" for all the cases which corresponds to 
the spectral resolution power of i? ~ 250 or ~ 30 
A at 8000 A for the 150 grism. The dispersion is 
2.8 A per pixel. Exposure time is typically 3600 
second. Unfortunately, these masks were not al- 
ways observed under the homogeneous condition; 
some exposures were obtained under cloudy /cirrus 
sky 

The data was reduced by using IRAF reduc- 
tion package. After bias subtraction and flat 
fielding, each slitlet spectrum was extracted, and 
then wavelength and flux-sensitivity calibrated. 
The flnal wavelength determination accuracy is 
a few A. More detailed description of the FO- 
CAS spectroscopy for SXDS sources is presented 
in Akiyama et al. (2005, in preparation). 

3.2. Results for OPEGs 

Table 3 summarizes the results for the OPEGs. 
At z' < 22, we observed 93 objects and the redshift 
of the 77 red galaxies (OPEG) were determined. 
65 arc in the 'general field' and 12 are in the 'high 
density' regions. At 22 < 2;' < 24, our complete- 
ness in determining redshift becomes low and the 
redshift of 22 galaxies in the 'general field' were 
obtained. Table 4 lists the results of the redshift 
measurement. 

Fig. 8 shows the spectra for the 77 red galaxies 
with z' < 22 in order of the obtained redshift. Ver- 
tical lines show the wavelength of Mgll A2800 (red 
dashed line in Fig.8), Fel A3581 (red), [Oil] A3727 
(blue), Mgl A3835 (red), Ca H and K (green), m 
(cyan), H7 (cyan), H/3 (cyan) lines and G band 
(green). For the most of the cases we see 4000 
A continuum break and Ca HK lines clearly and 
found that their spectra arc dominated by old stel- 
lar population. At the same time, [Oil] emission 
line with the observed equivalent width EWobs > 
A is detected in 29/77 objects. Only four galaxies 
have the rest-frame equivalent width larger than 
10 A and the strongest one has 16 A. Many of these 
lines are thus weak compared with the emission- 
line galaxies even in the rich clusters (Poggianti et 



al. 1999) but they are thus not evolving passively 
in very strict sense. The objects with [Oil] emis- 
sion are plotted in Fig. 13 below and we found 
that many of them are among the i' — z' bluer 
galaxies. There are also galaxies with weak but 
significant Balmer absorption lines (Fig.8). 

Fig. 9 shows the histogram of the obtained red- 
shift of the 77 red galaxies z' < 22 (solid line) as 
well as the in total of 99 red galaxies with z' < 24 
(dotted line). The redshifts were determined by 
the lines in each spectrum. The contribution of 
the 12 red galaxies in the high-density regions in 
Kodama et al. (2004) is also shown by the dashed 
line. We note that there is a strong concentration 
of the objects between z = 0.85 and z = 1.1 for 
galaxies with z' < 22. Indeed, 73 of the 77 objects 
with redshift (95%) is lie between z = 0.87 and 
1.12. The trend is also true for the fainter galax- 
ies although the fraction of the galaxies at z > 1.1 
is higher. Sampling galaxies in high density re- 
gion does not cause any difference (see below for 
the reality of the structure). Fig. 10 presents the 
same result in a different way; the redshift values 
are plotted against the z'-band magnitude. 

In Fig.ll, we plotted the observed i' — z' color 
of the objects along the redshift. At 0.8 < z < 1.1 
there is strong positive correlation between the 
i' — z' color and redshift. This is indeed expected 
for galaxies following the passive evolution. Wc 
plotted the passive-evolution models with Zf=2 
and 5 and L=5Lh< and O.IL, in the figure. Al- 
though the scatter exists, the data and the models 
agree quite well. From this figure, we strongly ar- 
gue that these galaxies with redshift indeed obey 
passive evolution and their red i' — z' color is due 
to the large 4000 A break of the old stellar popu- 
lation. Our color selection of OPEGs thus works 
very well. 

There is an obvious exception at z > 1.1; a 
galaxy at 2; ~ 1.4 has much redder color than the 
OPEG models. This is SXDS J021754.65-045027.6 
at z= 1.393. Since HK lines and the 4000A break 
are redder than our spectral range, we cannot say 
much about the nature of this galaxy or probably 
a reddened AGN. 

Of the 93 galaxies with 19 < z' < 22 selected 
by our color criteria, we could not determine the 
redshift of 16 objects. While it is true S/N ratio 
of spectra is not enough for some objects, which 
is due to poor sky condition, many of the spectra 
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of these objects arc rather featureless to securely 
determine their redshift. In Fig. 12, we plot one 
example of such red objects without any notable 
spectral feature to measure their redshift. 

We also obtained the spectra of non-OPEGs at 
0.85 < z < 1.1 in the same mask serendipitously. 
Fig. 13 shows the redhisft-z' — z' distribution of 
these galaxies as well as those of the OPEGs (the 
filled ciecles). Objects with EWob.s{OII) > 5 A 
are pitted with squares. Most of the bluer galax- 
ies have the [Oil] emission. This figure demon- 
strates how our selection criteria work to select out 
the passively-evolving galaxies. First, the OPEGs 
with [Oil] emission have the relatively blue colors 
above the criteria. It is then natural to consider 
that a fraction of OPEGs or their progenitors drop 
out from our criteria due to a little more amount 
of on-going star formation. On the other hand, 
some galaxies have bluer i' — z' color but no [Oil] 
emission. They may be the objects without sig- 
nificant on-going star-formation but with younger 
stellar population, such as the E + A galaxies seen 
in the intermediate-redshift clusters (e.g.. Dressier 
et al. 1999), or the younger OPEGs with forma- 
tion epoch less than z = 2. The fraction of the 
[OII]-detected objects is higher at ^ > 1 (19/36) 
than the lower redshift (10/41). This may be par- 
tially due to the selection effect since the colors of 
the OPRGs become redder as redshift increase and 
thus the rejection of the bluer galaxies are more 
sharp at the lower redshift. 

3.3. Galaxies in the high-density regions 

Fig. 14 shows the redshift distribution of the 
galaxies in the high-density regions of red galaxies 
in Kodama et al. (2004). The sample is dominated 
by the objects whose Ri'z' colors are consistent 
with the model galaxies at z ^ 1 with various ratio 
between old and young stellar population. Despite 
our speculation that there may be rich clusters of 
galaxies in these high- density regions, especially 
at CI, the redshift distribution does not show any 
strong single peak or concentration. CI and C2 
arc likely to be a projection of a few groups. CI 
has a rather discrete peak at z = 1.05. 

We also note that this result does not affect our 
conclusion in Kodama et al. (2004) since these 
high density regions of red z 1 galaxies were se- 
lected in order to get the statistical field subtrac- 
tion method work, and it does not really matter 



whether it is a single rich cluster at z 1 or a 
projection of a few groups at similar redshifts. 

4. NUMBER COUNTS 

As described in Sec. 2 and Sec. 3, we constructed 
a sample of OPEG candidates at z = 1 with 
Zf = 2 — 10 by the carefully determined color and 
magnitude criteria. Spectroscopic results, espe- 
cially for brighter galaxies, confirmed that the se- 
lection is highly efficient and at least 78% of the 
entire sample or 95% of those whose redshift were 
measured are indeed galaxies at 0.85 < z < 1.1 
whose spectra are dominated by old stellar pop- 
ulation. If we add a small number of objects at 
1.1 < z < 1.4, the fraction of OPEG is more than 
83% of the entire color-selected sample. Prom this 
fact, we may argue that our photometric sample is 
robust as it is dominated by OPEGs (we thus here- 
after refer the color-selected objects as 'OPEGs' 
instead of 'OPEG candidates' for convenience). 

In this section, we then show the magnitude- 
number relation of these OPEGs and their field- 
to-field variation at large scale. The luminosity 
function after taking the selected redshift range 
into account is then discussed in the next section. 

4.1. The average number counts 

Fig. 15 shows the number counts of the whole 
z'-band selected galaxies and as well as OPEGs. 
For OPEGs, we also plotted the counts after cor- 
recting the incompleteness (Fig. 6) evaluated for 
model galaxies with de Vaucouleurs light profile. 
For comparison, z'-band count obtained by Ca- 
pak et al. (2004) with the same instrument and 
filter is plotted. Although they agree with reason- 
able level, there is still slight (a few to several %) 
systematic difference between the counts, which 
corresponds to the w 0.05 mag zero-point differ- 
ence. This may be due to the field-to- field varia- 
tion or the different method in evaluating the total 
magnitude; Capak et al. (2004) used 3" -diameter 
aperture magnitude and aperture correction factor 
while we used the SExtractor MAG BEST values 
directly. 

While the z'-band number counts keep rising 
even at z' = 25, the coimts for OPEGs begin to 
flatten at z' ~ 22 and decrease toward the fainter 
magnitude even after the correction for the incom- 
pleteness. This is not due to the underestimation 
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of the incompleteness since, at z' ~ 24, OPEGs are 
in fact relatively more compact than the whole z'- 
selected ones and suffer less from detection incom- 
pleteness. The mode of FWHM values for OPEGs 
is 4.5 pix while that of the whole z'-selected ob- 
jects is 4.9 pix. Thus, while the number of the 
^'-selected sources which arc more diffuse than 
OPEGs still increase, the number of OPEGs de- 
crease along the magnitude. For the cases of pure 
passive-evolution models, the L, galaxy at z = 
is to be z' « 22 at 2 = 1. Therefore, if the se- 
lected red galaxies are dominated by OPEGs, the 
flattening at around z' — 22 is what wc expect for 
the case due to the shape of the galaxy luminosity 
function. 

4.2. Field-to-fleld variation at large scale 

We then investigated how the counts vary from 
fleld to field at large scale. In Fig. 16, we plot the 
number counts of OPEGs in each SXDS-C, -N, - 
S, -E, -W fields after the correction for the incom- 
pleteness. We found quite a large variation among 
them, especially in SXDS-S; at z' = 23, the high- 
est counts, of SXDS-N, is about twice as large as 
the lowest counts of SXDS-S. Note that the evalu- 
ated completeness at z' = 23 is 90-92% in SXDS-N 
and -S, respectively. The z'-band detection limit 
is also comparable, 25.7 and 25.5 for SXDS-N and 
-S, respectively. Therefore, we believe that the dif- 
ference of the number counts of OPEG candidates 
is not an artifact but the real nature of galaxies. 
We checked that the total counts of the z'-selcctcd 
galaxies have the much less scatter distribution 
among the fields, at most 5% (Table 2). 

To confirm that it is not due to the inho- 
mogeneous sample selection due to a systematic 
color difference among the fields, we also compare 
the R — z' color distribution of the galaxies with 
0.8 <i' - z' < 1.2 in Fig.l7. While the locations 
of the peaks of the color distribution of the five 
fields at aroimd i? — z « 1.85 agrees very well, a 
significant deficit of the galaxies in SXDS-S is seen 
at 1.6 < i? — z < 2.0, namely the color range of 
OPEGs. 

It is remarkable that the counts of OPEGs 
varies significantly even at this large scale. The 
each Suprime-Cam field of view has approximately 
34' X 27', which corresponds to 46 h'^Q Mpc x 36 
Mpc in comoving scale at z=l. The nominal 
comoving volume at 0.85 < z < 1.1 per field is 



1 X 10'^ Mpc^. The distribution of the mas- 
sive oldest galaxies appears inhomogeneous at 100 
Mpc scale at z = 1. 

In Fig. 16, not only the heights but also the 
shapes of the number counts show variation among 
the fields, in particular, the counts in SXDS-S 
show a rather flat slope at the faint end, while 
the counts in SXDS-N show a significant peak at 
around z' = 23 and a negative slope at the faint 
end. The flat slope in scon in SXDS-S may be 
caused by relatively larger contamination, since 
the overall number density of OPEG is the lowest 
in this field. The blue galaxies, which are abun- 
dant, tend to be scattered into our color cut due 
to their photometric errors. Also, the background 
galaxies that have similar colors as those of the 
2=1 OPEGs inevitably come into our sample es- 
pecially at the faint end. We repeated the similar 
procedure in Sec. 2. 2. 4 in the SXDS-S region and 
found that the contamination due to the random 
photometric error is larger in this field compared 
to other fields, but is still only 2% for the objects 
at z' = 24-25. As for the contamination by the 
background galaxies, since our spectroscopic sam- 
ple is nearly complete only down to z' ~ 22, the 
true redshift distribution of the sample at the faint 
end is unknown. If the background contamination 
in the SXDS-S field is as large as « 50% at z' = 24- 
25, the slope of the faint end would become nega- 
tive and all the fields would have consistent slopes 
after subtracting the background. 

We note that Kodama et al. (2004) reported a 
much stronger deficit (a factor of 2-3 from z' — 22 
to 24) of red faint galaxies at z ~ 1 in the same 
SXDS field, even though they applied a similar 
colour selection to ours based on Riz' . As shown 
in our Fig. 14, the averaged number counts of the 
OPEG candidates over the entire SXDS field does 
not show such a prominent deficit towards the 
faint-end, although the slope is clearly negative. 
It should be reminded, however, that Kodama 
et al.'s (2004) sample of red galaxies are differ- 
ent from that in this paper in two ways: First 
of all, their sample were constructed in the se- 
lected high density regions of red z ~ 1 galax- 
ies, rather than the entire SXDS field used here. 
Secondly, they have applied a statistical subtrac- 
tion of foreground/background contamination us- 
ing the mean density of such red galaxies outside 
the high density regions. To discuss more on this 
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interesting issue on the deficit of faint red galaxies 
(which is often referred to as "down-sizing"), we 
badly need more accurate photometric redshift in- 
formation on these faint galaxies by adding deep 
NIR imaging which will be provided in the near 
future through the UKIDSS collaborations. 



a 



o 









— 






































1 

< 


1 












































































SXDS 


Centt 


































SXDS 
SXDS 


Nq 

So 


rtt 

utl 


1 








* 
























SXDS 


East 
























































SXDS 
SXDS 


-Wi; 
Al 


st 


/e 



















































































19 20 21 22 23 24 25 26 



z [AB. Magbesl] 



Fig. 6. — Completeness of the detection of OPEG 
candidates on the z' band images in SXDS-C, - 
N, -S, -E, -W fields evaluated by using the model 
galaxies and the observed background-noise im- 
ages. De Vaucouleurs profile is assumed for the 
model galaxies. 
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5. DISCUSSIONS 



5.1. Luminosity Function of OPEGs at z = 

1 

Provided the redshift selection function based 
on the spectroscopic sample, we discuss the lumi- 
nosity function of the OPEGs and their spatial 
density. We here use the simplest procedure to 
obtain the rest-frame i?-band luminosity function 
(LF) of OPEGs at z = 0.9-1.1 based on the result 
that the redshift distribution of the z' < 22 spec- 
troscopic sample show strong peak at the redshift 
range. 

With minimum assumption, we locate all the 
selected OPEGs at z = 1 and adopt the single 
value of the fraction of galaxies between z=0.9 and 
1.1 from the spectroscopic results. The difference 
in distant modulus between z — 0.9 and 1.0, and 
that between z = 1.0 and 1.1 is 0.28 and 0.26 mag, 
respectively. The analysis yields magnitude uncer- 
tainty of these values, which is still smaller than 
the bin width of 0.5 mag adopted here. For the 
spectroscopic sample of 93 z' < 22 OPEG candi- 
dates, 59 are at z = 0.9 — 1.1, which yields 63.4%. 
We simply calculate absolute magnitude a,t z — 1.0 
for all the OPEG candidates and multiply the fac- 
tor 0.64 to the density values at all the magnitude 
bin. Fig. 18 shows the results. We first limit the 
magnitude range to Mb = —21 which corresponds 
to z' = 22.2 at z = 1.0 and is close to the depth 
of the spectroscopic sample. The red circles are 
the LF obtained from the raw counts and the blue 
ones show that after the incompleteness correction 
using the whole-field data. We discuss the differ- 
ence from the LF of the local early-type galaxies 
in the next subsection. 

We then extend the analysis to the true limit 
of the photometric sample, z' = 25 as shown in 
Fig. 19 assuming that the redshift selection func- 
tion and the contamination fraction are the same 
as the brighter objects above our spectroscopic 
limit. We note that the true fraction of con- 
tamination may increase towards fainter magni- 
tudes due to the background galaxies whose col- 
ors match with our criteria. For example, the old 
red galaxies at z > 1.1 with some right amount 
of star-formation could have the similar color as 
the OPEGS at z = 0.9-1.1. Thick line is the 
best-fit Schechter function for the incompleteness- 
corrected data. Thin lines arc the same but for 
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Fig. 7. — The sky distribution of the targets in 
the spectroscopic observations. The circles are the 
OPEG candidates and the squares are other ob- 
jects. Filled symbols are those for which redshift 
were obtained. Large circles show the high-density 
regions of red galaxies studied in Kodama et al. 
(2004). 



12 



40 - 

.1 






^A* ' ■ 


; ; z= 


>1 


20 



6000 BOOO 

joai'Boa.BO-;045e3i e ^|^^^-4o 

6000 BOOO 


6000 

r jo2ie34.Bb- 

6000 


8000 

isoHo.'a . 
eooD 


6000 

L; Jp3 1833,02-: 050 
: ■ z= 0.B73 ■ ■ ■ ■ 

6000 


8000 
8000 


40 

»« 


J021B04.85- 


0a9913B k-40 


J0Eia09.£5- 


43449 'a jli4 


L; J021824.26-:044 
: ; 7,= ( H!JO : : : : 


1 - 


f ao ■ 


J0H1947.0E- 
6000 


BOOO 
BOOO 


6000 
J0ai939.74- 

6000 


800D 
BOOO 


L ; JOE 1 939,9 t'-D5QB0E.l^^^^^^^^ 

fio 111 inimi 


40 - 
so - 



joai656.3a- 

1= 0.895 : 




jossoie. ■ 

z= 0,913 : 







aooo aooo eooo aooo (iiioii iiooo 

Observed Wavelength [A] Observed Wavelength [S] Observed Wavelength [S] 
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Fig. 9. — The redshift distribution of the spectro- 
scopic sample of OPEGs. 
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each SXDS-C, -N, -S, -E, and -W fields; density 
variation is seen as expected from the surface num- 
ber counts. The resultant Schechter function pa- 
rameter values are summarized in Tabic 5. 

For the check, we also calculated the LF from 
the sample selected by the 3" -diameter-aperture 
colors in order to see how large is the uncertainty 
of the color selection. Indeed, the i' — z' color 
in 3"-aperture is ~ 0.05 mag bluer than the 2"- 
aperture. The results are shown in Fig. 20 and Ta- 
ble 5. 3" -aperture color selection yields slightly 
but systematically smaller number density and 
steeper faint-end slope although the difference be- 
tween the Schechter parameters of the two cases is 
small and just marginally significant. The differ- 
ence is due to the fact that the galaxy colors move 
systematically bluer and a fraction escaped from 
the model criteria (Fig. 4). 

5.2. Comparison with local ecirly-type 
galaxies 

The colors of the local counterparts of 2 = 1 
OPEGs can be inferred by extrapolating the pas- 
sive models to z = 0. Although secondary bursts 
due to minor merger events may change their color 
in short time; sc;ale. it does not afFcicts the spec- 
tra of the entire galaxy after a few times 10* yrs 
(see Sec. 2. 2.1). Thus the descendants of 2; = 1 
OPEGs should be among the local galaxies which 
form the red color-magnitude sequence expected 
for old (> 10 Gyr) passive galaxies. 

We here compare the number densities of 2: = 1 
OPEGs with the red galaxies detected in Slone 
Digital Sky Survey (SDSS). To avoid any system- 
atic zero-point error, we consider the color range 
from the observed red sequence of the local galax- 
ies instead of directly extrapolating the models to 
define the OPEGs at z = 1. The color-selected 
r'-band LF of the SDSS red-sequence galaxies is 
kindly provided by Osamu Nakamura (in private 
communication). We referred the red sequence 
(zero points) of the morphologically-selected E/SO 
galaxies in Fukugita et al. (2004) and adopted the 
color cut at the blue side A(g' — r')=— 0.1. The r'- 
band LF was calculated as in the similar manner 
in Nakamura et al. (2003) and is shown in Fig. 18 
by the thick dashed line; we assume a single rep- 
resentative color, B — r'=1.2, to shift it for the 
B-band. 



It is also interesting to compare z = 1 OPEGs 
and the morphologically-selected E/SO galaxies 
at 2; = themselves since most of the local 
early-type galaxies have no or only a very small 
amount of young stellar population and their spec- 
tra can be well explained by the passive evolu- 
tion models with old age (e.g., Kodama and Ari- 
moto 1997). We compare our result with the LF 
of the morphologically-selected early- type galaxies 
by Nakamura et al. (2003), which is also plotted 
in Fig. 18 by the solid thick line. The magnitude 
shift, B — r'=1.2, is also adopted. The number 
density of the color-selected red galaxies in SDSS 
is in fact larger than that of the morphologically- 
selected ones. The difference may be due to the 
presence of the morphologically later-type galaxies 
with similarly red colors. 

It is clear in Fig. 18 that the number density 
of the OPEGs more Imninous than Mb = —21 
is much larger than those of the local 'early-type' 
galaxies. 

We plotted in Fig. 18 the evolved luminosity 

function; pure luminosity evolution for passive 
models with Zf = 2 and Zf = 5 are plotted with 
the cyan and the green lines, respectively. We 
used the representative magnitude shift values, 
AMb — —1.5 and —1.1 for the cases of 2/ = 2 
and 5, respectively. The prediction with Zf > 3 
is similar {AMb < 0.1 mag) with the case of 
Zf = 5. Since our sample of OPEGs is selected 
by their colors based on the similar passive evolu- 
tion models, we can directly compare the number 
densities of the observed results with these pure 
passive models in a very consistent manner. The 
number density of OPEGs at z '--^ 1 is « 65, 50, 
and 30% of the predicted ones of the evolved lo- 
cal color-selected red galaxies and w 90, 85 and 
85% of the local morphologically-selected galaxies 
at the evolved absolute magnitude Mb — — 1, 
M*, M, -I- 1, respectively, for the case of 2/ > 3. 
If the local early-type galaxies are as young as 
Zf = 2, OPEGs occupy only « 20-40% and « 40- 
80% of the progenitors of the local color-selected 
and morphologically-selected early-type galaxies, 
respectively, between Mb — M, — 1 and -|- 1. 
As it is not realistic that all the local early-type 
galaxies are as young as 2/ = 2, as we see at least 
some evolved galaxies at 2; > 2 (Daddi et al. 2005), 
it is fair to say that more than half, but not all, of 
the present-day typical (M* - 1 < M < M^ + 1) 
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carly-typc galaxies had already been formed into 
luminous quiescent galaxies at z = 1. 

The results also mean that the less than half, 
but non-negligible fraction of the present-day 
early- type galaxies have been evolved since ^ ~ 1. 
There are at least three possibilities to explain 
the result; more star-formation at z ~ 1, the 
younger age (they are quiescent but formed at 
1 < z < 2), and the evolution via merging. Some 
different star-formation history to intermediate 
redshift may be allowed so that they are observed 
as OPEGs at 2 = if there are enough long time 
(~ 5 Gyr) after the star- formation ceased (e.g., 
Bower, Kodama, Terlevich 1999). That is also 
true for the case that galaxies are younger than 
z = 2, i.e., the major star- formation finished be- 
tween z = 1 and 2. They are too blue to be 
classified as the OPEGs defined here but again 
there is enough time to be the OPEGs a,t z = 0. 
In either case, the progenitors can be among the 
different category of the objects (e.g., late- type 
morphology). There is also a possibility that a 
fraction of the early-type galaxies at 2 = have 
experienced a major merger event as other recent 
observational results also suggest (e.g., Bell et al. 
2005; van Dokkum et al. 2005). Even for the case 
of the 'dry merger' (Bell et al. 2005), the resul- 
tant luminosity shift of factor ^ 2 can explain the 
difference between our z = 1 OPEG LF and the 
LP of 2; = morphologically-selected early-type 
galaxies. 

On the other hand, our results show that stars 
in more than half of the OPEGs at z = 1 were 
formed by z ~ 2. It can be tested by observing the 
OPEG LF between z = 1 and 2 whether such pro- 
cess as 'dry merger" dominantly occurred or not; if 
we see much less number density of OPEGs than 
predicted by the passive evolution models, that 
is to be evidence that galaxy assembly occurred 
later than their star formation for the early-type 
galaxies in general field. 

5.3. Comparison with other results at in- 
termediate redshift 

We also compare our results with the number 
densities of old or early- type galaxies at intermedi- 
ate redshift selected with various different criteria. 

There were many results about the number den- 
sity of color- or morphologically-selected early- 



type galaxies at z ^ 1 (Lilly et al. 1996; Kauff- 
mann. Chariot, & White 1996; Zepf 1997; Totani 
& Yoshii 1998; Prancheschini et al. 1998; Berger 
et al. 1999; Sehade et al. 1999; Menanteau et al. 
1999; Treu & Stiavelh 1999; Thompson et al. 1999; 
Yan et al. 2000; Kajisawa and Yamada 2001; Im et 
al. 2002; Daddi et al. 2002: Cimatti et al. 2002a, 
b, c; Roche et al. 2002; Miyazaki et al. 2003; 
Pozzetti et al. 2003; Glazebrook et al. 2004; Bell 
et al. 2004; Saracco et al. 2005; Treu et al. 2005). 
Treu (2003) present an excellent review for these 
previous works. The samples of galaxies in these 
previous study, however, are limited to relatively 
small area of at most several hundred arcmin^ or 
to the small numbers of at most a few hundred. 

Miyazaki et al. (2003) obtained the luminos- 
ity functions of their sample of Extremely Red 
Objects. They first selected EROs by the crite- 
ria R — Kab > 3.3 and then classified them to 
'old galaxies' (OGs) and the 'dusty galaxies' by 
the analysis of the multi-band SED data. They 
found that the luminosity function of the OGs 
are consistent with the passive evolution predic- 
tion from the LF of the z = morphologically- 
selected E/SO galaxies by Marzke et al. (1998). In 
fact, our OPEG LF also coincides with the passive- 
evolution prediction using the local LF of Marzke 
et al. (1998) but Marzke et al. (1998) gives lower 
number density compared with Nakamura et al. 
(2003). In this sense, our results seem consistent 
with the LF by Miyazaki et al. (2003) but at the 
same time the conclusion is different and that the 
LF of early-type galaxies do show evolution if we 
refer the local LF by Nakamura et al. (2003) based 
on the homogeneous large sample using the SDSS 
database. 

Bell et al. (2003) reported their results of the 
COMBO- 17 survey of 0.78 deg^. The multi-band 
photometric data allow them to evaluate photo- 
metric redshift of the galaxies and to make the 
robust spectral classification in the volume compa- 
rable (60%) with SXDS. The limiting magnitude 
of the sample is ~ 22 AB magnitude at 9140 A , 
which is roughly equal to an OPEG with L» lumi- 
nosity at z = and 3 magnitude shallower than 
the sample in this paper. They found the strong 
color bimodality among their sample and the mean 
color of the red galaxy sequence evolves with red- 
shift in a way that is consistent with the aging 
of old stellar population. At the same time, their 
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B-band luminosity function of the rcd-scqucnce 
galaxies shows very strong evolution toward z = 1. 
The density normalization 0* at 2; ~ 1 is almost 
15% of the value at z = 0. In our data, however, 
no such strong evolution in the density normal- 
ization is seen. At Mb = —21, the density of 
the red-sequence galaxies at z = 0.95 in Bell et 
al. (2004) is only « 50% of the OPEGs at z = 1 
in this paper. The difference may not be due to 
the difference of the color-selection criteria. While 
Bell et al. first fitted the color-magnitude sequence 
(based on the photometric redshifts) and then take 
the galaxies down to 0.25 mag bhicward in rest 
U — V color, our criteria sample the galaxies to- 
ward « 0.15 mag bluer R — z' (close to the rest 
U — B) color from the peak (Fig. 4). For the typi- 
cal galaxy SED, A{U — V)= 0.25 corresponds to 
A(f7 — B) w 0.15. Sample completeness or field- 
to-ficld variation should be examined further to 
understand the difference. 

Chen et al. (2003) studied the number den- 
sity of color-selected early-type galaxies using the 
if-band selected sample of the ~ 0.35 deg^ area. 
Their detection limit is H = 20.8 (5(j). Since their 
color classification of red galaxies include not only 
E/SO galaxies but also Sab galaxies, it is difficult 
to directly compare their result with ours. They 
obtained a result that the number density of these 
early- type galaxies decreases ~ 40% since 2=1. 

Among these studies, the current sample of 
SXDS z — 1 color-selected early-type galaxies, or 
OPEGs, is the largest in both the sample size (the 
number of galaxies) and the survey area, which 
must considerably increase the statistical accuracy 
of the results. Moreover, we made a pilot spectro- 
scopic observation to calibrate our 'photometric' 
sample, although the number of the galaxies with 
redshift information is still limited. 

6. Summary 

We have constructed a large sample of the » 

4000 color-selected old passively-evolving galaxies 
(OPEGs) at z = 0.85 — 1.1 using the notably deep 
and wide-field imaging data of the Subaru/XMM- 
Ncwton Deep Survey. Although our sample selec- 
tion was first made by the simple color-cut based 
on the galaxy-evolution models but the obtained 
color and the magnitude distributions as well as 
the results of the spectroscopic observations of 



~ 80 galaxies strongly supports that the sample is 
indeed dominated by the OPEGs at z = 0.85 — 1.1. 
The number-density of the OPEGs still show 10- 
30% variations over the 0.5 deg scale and the whole 
average density is « 40 — 50% of those of the local 
color-selected red galaxies and 60-85% of the local 
morphologically-selected galaxies at Mb ~ Af*. In 
the near future, deep near-infrared data from the 
UKIDSS project as well as the mid-infrared data 
from the SWIRE survey provides further complete 
SED coverage and the whole evolutionary feature 
of OPEGs can be studied over large redshift-range. 

We thank the staff of the Subaru Telescope for 
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entific research of the Ministry of Education, Cul- 
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Table 1 

Summary of the SXDS Catalog 



Field 


Aclcan 


FWHM 


Tcxp(J?) 


Tcxp(^ ) 


Tcxp(2') 


Rlim 


./ 

*lim 


/ 

^lim 




arcmin^ 


arcsec 


sec 


sec 


sec 


mag 


mag 


mag 



SXDS Center 


960.1 


0.80 


14880 


38820 


13020 


27.0 


27.1 


26.0 


SXDS North 


938.2 


0.84 


13920 


26405 


15122 


27.2 


27.1 


26.1 


SXDS South 


955.8 


0.82 


13920 


18540 


11040 


27.1 


26.9 


25.8 


SXDS East 


905.3 


0.82 


13920 


22060 


15990 


27.0 


26.9 


25.9 


SXDS West 


904.4 


0.82 


13920 


35883 


18660 


27.0 


27.0 


26.1 




Fig. 11. — i' — z' color and redshift distribution of 
the OPEG candidates with spectroscopic redshift. 
Passive evolution models for zp = 5 (upper) and 
2 (lower) and 5L* (thick) and O.IL* (thin). 



20 



Table 2 

The Number of the OPEG Candidates 



Field 


Aused 


vl{z' < 25) 


noPBG(a' < 25) 




arcmin^ 






SXDS Center 


832.3 


39166 


987 


SXDS North 


811.1 


39670 


981 


SXDS South 


830.2 


38138 


653 


SXDS East 


781.2 


37201 


1016 


SXDS West 


778.3 


35629 


843 


SXDS All 


3690.0 


175601 


4118 



Table 3 

Summary of the Spectroscopic Observations 





field 
n(with z) 


n(no z) 


HDR 
n(with z) 


n(uo z) 


19 < / < 22 


65 


12 


12 


4 


22<z' < 24 


22 


58 





2 


total 


87 


70 


12 


6 
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Table 4 
Redshifts of z' < 22 OPEGs 
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Table 4 — Continued 



Name Zbcst i' - z' R-z' Redshift" EWobs{OII) 

(mag) (mag) (mag) A 



J021736, 


.51- 


-051035, 


,9 


21 


.65 





.93 


1 


.92 


0.987 


5 
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,9 


20 


.95 
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1 


.94 


1.002 





J021718, 
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.4 


21 
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1 
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1.009 
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J021708, 


.79- 
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.7 


20 
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1 


.03 


1 
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1.015 
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,6 


21 
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1 
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1 
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21 
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1 
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,2 


21 
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1 
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21 
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1 


.84 


1.033: 
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20 
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20 


.47 
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1 
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20 
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1 
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21 
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.01 
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1 


.78 
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.1 


21 
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1 
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1.060 
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,9 


21 


.53 


1 
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1 
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J021734, 
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.4 


21 


.42 
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2 


.02 


1.063 
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.4 
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.39 
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1 


.86 
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21 


.59 
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.00 
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1.074 
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.93 


1.075 





J021746. 


.30- 


-045211, 


.1 


21 
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.43 
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.02 
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.94 
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,2 
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.95 


1 


.81 
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,8 


21, 
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1, 
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1 


.95 
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.84 
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J021650, 
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,8 


21 


.97 
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1 


.81 
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J021800, 


.18- 


-050319, 


,6 


21 
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20 
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2 
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21 
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.95 


1 


.88 


1.096 
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J021946, 


.82- 


-050633, 
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21 
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.11 


1 


.93 


1.120 
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,8 


21 


.69 


1 


.07 


1 
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1.186 
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-044757, 
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21 


.40 
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.04 


2 
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,4 
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1 
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1 
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24 



(a) Rcdshift values with ':' have uncertainty 
due to the less prominent spectral features. 
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Table 5 

schechter fitting parameters for the b-band luminosity function of opegs 







err. 


4>* 


err. 


a 


err. 




AB mag. 




xlO"^ Mpc"^ 








All (single-z method) 


-21.38 


0.10 


1.46 


0.14 


-0.67 


0.07 


Center 


-21.34 


0.21 


1.52 


0.32 


-0.72 


0.15 


North 


-21.29 


0.21 


1.74 


0.30 


-0.60 


0.15 


South 


-21.33 


0.27 


1.02 


0.24 


-0.63 


0.18 


East 


-21.60 


0.21 


1.46 


0.28 


-0.76 


0.14 


West 


-21.39 


0.22 


1.46 


0.32 


-0.69 


0.16 


All (3" aperture) 


-21.25 


0.11 


1.38 


0.12 


-0.54 


0.09 
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Fig. 12. — An example of the objects with feature- 
less red continuum. 



Fig. 14. — Redshift distribution of the galaxies in 
the 'high density region' no.l and 2 in Kodama et 
al. (2004). 
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Fig. 13. — Color distribution of the all the galax- 
ies at 0.8 < z < 1.15. The filled circles arc the 
OPEGs. The open squares show the objects which 
have the [Oil] emission lines with the observed 
equivalent width EWobs{OU)> 5 A. 
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Fig. 15. — Number counts of the z'-selected ob- 
jects (top) and the OPEG candidates (bottom). 
The values after incompleteness correction for the 
OPEG candidates are shown by the open circles. 
The z'-band number counts obtained by Capak et 
al. (2003) using slightly different method is also 
plotted by the open squares for comparison. 
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Fig. 16.— Number counts for OPEGs in SXDS-C, 
-N, -S, -E, -W fields. 
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Fig. 18.— z = 1 OPEG luminosity func- 
tion at Mb < —21 is shown. The red and 
blue circles show those before and after the cor- 
rection for the detection incompleteness. For 
comparison, the local luminosity functions for 
the SDSS morphologically-selected E/SO galaxies 
(solid black line, Nakamura et al. 2003) and for the 
color-selected red galaxies (dashed black line), as 
well as those at z = 1 as expected by pure lumi- 
nosity evolution with the same passive-evolution 
models with Zi?=2 (blue lines) and 5 (green lines) 
are plotted. The vertical lines show the Af* values 
for the z = E/SO galaxies. 
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Fig. 17. — R-z' color distribution for the OPEGs 
in SXDS-C, -N, -S, -E, -W fields. 
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Fig. 19. — z = 1 OPEG luminosity function. 
The open (red in the on-Une edition) and the 
filled (blue) circles show those before and after the 
correction for the detection incompleteness. The 
Schechter function fit for those obtained in each 
SXDS-C, -N, -S, -E, -W fields (thin lines) are plot- 
ted together. 
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Fig. 20. — z = 1 OPEG luminosity function (smae 
as in Fig. 19). The Schechter function fit for the 
3"-aperture case (thin line) is plotted together. 
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